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A mapping method is developed to investigate the problem of determination 
and control of heat-deposition patterns on the plates of a tokamak divertor. 
The deposition pattern is largely determined by the magnetic field lines, which 
are mathematically equivalent to the trajectories of a single-degree-of-freedom 
time-dependent Hamiltonian system. Maps are natural tools to study the 
generic features of such systems. The general theory of maps is presented, and 
methods for incorporating various features of the magnetic field and particle 
motion in divertor tokamaks are given. Features of the magnetic field include 
the profile of the rotational transform, single- versus double-null divertor, 
reverse map, the effects of naturally occurring low M and N, and externally 
imposed high-M, high-N perturbations. Particle motion includes radial 
diffusion, pitch angle and energy scattering, and the electric sheath at the plate. 
The method is illustrated by calculating the stochastic broadening in a single- 
null divertor tokamak. Maps provide an efficient, economic and elegant method 
to study the problem of motion of plasma particles in the stochastic scrape-off 
layer. 


1. Introduction 


Divertors are a common feature of present-day large tokamaks. The heat 
fluxes striking the collector plates may be as high as 15-20 MW m 1. Issues 
related to the heat flux deposition are of considerable importance for the design 
of future tokamaks. Therefore the prediction and preferably the control of the 
width of the scrape-off layer of a divertor is a central issue in the design of 
ignition tokamaks such as the international thermonuclear experimental 
reactor ITER (Post et al. 1991). Divertor simulations generally assume that the 
plasma is bounded by a sharp separatrix, which is correct only if the toroidal 
symmetry is perfect. Toroidal asymmetries in the magnetic field cause the 
toroidal magnetic flux inside the last confining magnetic surface to be smaller 
than the flux inside the unperturbed separatrix, and create a layer of open 
magnetic field lines that lies between the last confining magnetic surface and 
the unperturbed separatrix. This layer is called the stochastic layer. It is 
important to know the effects of naturally occurring perturbations, such as 
toroidal ripple or magnetohydrodynamic (MHD) modes, on the design of 
divertors. The introduction of magnetic perturbations may allow an opti- 
mization of the width of the divertor scrape-off layer and control of the heat 
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deposition on the divertor plates. Thus the issue of heat deposition on divertor 
plates is of critical importance in tokamak physics. 

We have developed a new approach to investigate this problem that makes 
use of what we call tokamak divertor maps. This approach is based on the 
observation that the magnetic field trajectories in a divertor tokamak are 
mathematically equivalent to a single-degree-of-freedom, time-dependent 
Hamiltonian system, and that the basic features of motion near a separatrix 
broadened by asymmetric perturbations are generic for such Hamiltonian and 
near-Hamiltonian systems. In this approach, the motion of field lines, including 
the effects of asymmetric perturbations, is represented by area-preserving maps 
and their continuous analogues, which we call tokamak divertor maps. The key 
results obtained using this method to calculate stochastic broadening in a 
single-null divertor were reported by Punjabi, Verma & Boozer (1992). In this 
paper we give a thorough exposition of the method, with a detailed discussion 
of how it is applied to the problem of determining the pattern of the strike 
points of particle trajectories on the divertor plates of a single-null divertor 
tokamak. 

Pomphrey & Rieman (1992) have made related calculations using wires to 
model divertor configuration instead of a map. B. J. Braams & M. D. Teytel 
(1993 personal communication) have employed the grid-free random-walk 
method to study this problem. La Haye’s (1991) work uses a differential- 
equation approach with a realistic magnetic geometry for the DIII-D device. 
Our mapping method, on the other hand, is much more computationally 
efficient and suitable for looking at detailed generic issues. 

The remainder of this paper is organized in six sections. In §2 we state the 
purpose for developing this method, and describe the desirable features that we 
should like it to possess. In §3 we introduce the idea of Hamiltonian mechanics 
of a broadened separatrix. In §4 we introduce the simplest of the tokamak 
divertor maps, called the simple map. We also present some useful extensions 
of the simple map. These include a map for the double-null divertor tokamak, 
a continuous analogue of the simple map, the backward map, a map with high 
adjustable shear, and maps that include naturally occurring low poloidal and 
toroidal mode numbers (M and N) as well externally imposed high-M, high-N 
perturbations. We also give two methods for calculating the deposition pattern 
of field lines on the divertor plates. In §5 we give the general theory of maps, 
which forms the theoretical foundation for the tokamak divertor maps. In this 
section, we also discuss the generating functions for area-preserving maps, and 
the explicit, reversed and continuous maps. In §6 we discuss a Monte Carlo 
procedure for calculating the heat deposition on divertor plates. This procedure 
includes the effects of spatial diffusion, velocity-space scattering, energy 
scattering and the electric sheath at the divertor plate on the particle 
trajectories. In §7 we show in detail how the simple map is applied to a single- 
null divertor tokamak to calculate stochastic broadening of the scrape-off layer. 
We also summarize and discuss our results in this section. 


2. Purpose of the method 


Our purpose is to calculate the effects of magnetic asymmetries on the heat- 
deposition pattern on divertor plates, and to determine if the heat deposition 
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can be broadened. We have developed a method called the method of maps to 
investigate this problem. This new method correctly addresses the issues of 
magnetic and particle effects. Magnetic effects involve the topology of the 
toroidal magnetic configuration, the width of the stochastic scrape-off layer, 
and the rotational-transform profiles. Particle effects include motion along the 
field lines, radial diffusion, and pitch angle and energy scattering. Since the 
mean free path of electrons in the stochastic scrape-off layer is comparable to 
the magnetic correlation length, a fluid treatment of this problem is not 
satisfactory. The method of maps is efficient, economical and flexible, and the 
physics of the method is transparent and clear. 

In the simulations, which assume a sharp separatrix, the width of the scrape- 
off layer is determined by the assumed diffusion coefficient, which is usually of 
order 1 m? s™!. The field lines of the scrape-off layer strike the divertor plates in 
a stripe with a width that is determined by axisymmetric magnetics. For 
asymmetric perturbations the width of the stochastic layer near the X point is 
approximately proportional to the square root of the toroidally asymmetric 
perturbation amplitude (Boozer 1983). If the magnetic perturbations are 
sufficiently large, the scrape-off layer is dominated by the stochastic layer. The 
field lines of the stochastic layer strike the divertor plate in a stripe that is much 
narrower (by a factor of about six) than the stochastic region of magnetic flux 
that lies between the last confining magnetic surface and the ideal separatrix. 
One of the important quantities used in this paper is the width of the stochastic 
scrape-off layer. This width is measured as the distance from the X-point to the 
last confining magnetic surface, where the distance is measured along the line 
from the O-point to the X-point. The greater width of the scrape-off layer in 
comparison to the divertor stripe has important implications for the design of 
divertors. 


3. Hamiltonian mechanics of a broadened separatrix 


Any globally divergence-free field in 3-space is mathematically equivalent to 
a single-degree-of-freedom Hamiltonian system (Whittaker 1937; Cary & 
Littlejohn 1983; Lichtenberg & Liberman 1983). For example, the motion of 
magnetic field lines in a tokamak can be described in terms of a poloidal angle 
0 as the canonical co-ordinate, toroidal flux / as the canonical momentum, the 
negative of poloidal flux, namely y = y(w,6,¢), as the Hamiltonian, and a 
toroidal angle ¢ as the canonical time. The phase space of this Hamiltonian 
system can be extended by considering the negative of the Hamiltonian, 
namely —y, and ¢ as conjugate variables. The phase space is then three- 
dimensional, and phase flow is parametrized by canonical time. For closed 
magnetic confinement systems, invariant surfaces are nested tori in 3-space. 
These Hamiltonian systems have the same topological structure as the simple 
map. The basic features of the motion of these systems near a separatrix are 
generic. Therefore the Hamiltonian mechanics of the broadened separatrix is 
the natural approach to investigate the motion in the scrape-off layer of a 
divertor tokamak. This approach can be made efficient and economic if we can 
construct area-preserving discrete maps to represent the Hamiltonian systems. 
The continuous analogue of the maps should be solutions of the single-degree- 
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Ficure 1. For caption see facing page. 
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of-freedom time-dependent Hamiltonian equations of motion. For our purpose, 
we want the maps and their continuous analogues to have toroidal geometry. 

First, we present an example of such maps, called the simple tokamak 
divertor map. 


4. Simple tokamak divertor map 
The simplest divertor map has a single positive parameter k, 


İni =2,—ky,(1—y,), (1) 
Ynrı = Yn + ktn (2) 


an O-point at x = 0, y = 0, and an X-point at x = 0, y = 1. The map provides 
a good representation of the magnetic field lines of a generic single-null 
tokamak divertor (figure 1a). Toroidal asymmetries can be simulated using k or 
by including additional terms. Here, toroidal asymmetries are simulated using 
k, but the results are not qualitatively changed by other representations. The 
smaller k is, the closer the last confining magnetic surface passes to the X-point. 
For k = 0:6 the closest approach is x = 0, y = 0-997 (figures La, b); for k = 09 
the closest approach is x = 0, y = 0:854 (figure 1c). Each iteration of the map 
corresponds to a toroidal angular advance ¢,. Near the O-point the rotational 
transform of the field lines per iteration, 1 = 27/q,, is 


sin (27/q,) = (1-34). (3) 


One obtains a transform of roughly unity near the axis if 10 iterations of the 
map constitute a single toroidal circuit, £, = 27/10. 


4.1. Extensions of simple map 


The simple map can be extended to represent various magnetic features of 
divertor tokamaks. These include double-null divertors, realistic transform (e.g. 
high shear) and more physical perturbations, etc. 

The magnetic topology of a double-null divertor tokamak can be adequately 
modelled by the map 


xı = X%y—ky,(1 =Y), (4) 
Yı = Yo + kay. (5) 


The rotational transform of divertor tokamaks can be reproduced if we consider 
Nmap © 10 periods of the simple map to be equivalent to one toroidal circuit of 
a field line in a tokamak. To plot the intersection of trajectories with the 
divertor plates, we require the location of trajectories at arbitrary toroidal 
locations, not just at the principal planes of the discrete map. This can be 
accomplished by constructing the continuous analogue of the simple map. 


FIGURE 1. (a) The simple map when the parameter k = 0-6. The figure shows the separatrix, 
the plane y = 1 and the divertor plates. (b) Magnified view of the stochastic scrape-off layer 
near the X-point when k = 0-6. This plot was done by iteration of the initial condition 
£o = 0, Ya = 09971. (c) Simple map when k = 09. Now the last good surface passes near 
x = 0, y = 0-854. The stochastic layer here is plotted from iterations with initial condition 
Ly = 0, Ya = 085492. 
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Letting ¢ = 6/f,,0<¢< 1, the (area-preserving) continuous analogue of the 
simple map is 
Aİ) = 2,—kyn(1—Yn) $, (6) 
Yİ) = Yn tkd). (7) 
Note that these equations are the solutions of Hamilton’s equations of motion 
with 
H(x,y) = — kx(x— xo) — kyyo(1 — yo), (8) 


where x is the canonical co-ordinate, y is the canonical momentum and ¢ is the 
canonical time. The continuous map is essential for finding the toroidal 
bunching of the location of field line strike points, as well for simplifying the 
Monte Carlo calculations for the heat-deposition pattern on the divertor plates. 

Particles that have a negative parallel velocity follow the map backwards. 
The reverse of the simple map is 


Yo = Yı — kz, (9) 
Ly = Tı tkyoll— yo). (10) 


The simple map has low shear away from the separatrix. By adding terms to 
the simple map, the shear may be increased. 


4.2. Higher shear map 
An explicit map that has an adjustable shear s is 


T = X— ky [(1—yo) (1+ sy) + 527), (11) 
y= Yo + ka, [i ts(x1* yol). (12) 


Tokamak shear corresponds to negative s. For k = 4 the ratio of ı at the X-point 
to 4 at the surface that goes through (x = 0, y = 0:8) is 1:35 for s = 0 but 211 
for s = — 0'8. Near the X-point at (x = 0,y = 1) the map has the same form as 
the simple map, but with k—> k(1 + s). Thus it is possible to construct maps with 
variable shear. 


4.3. Perturbations 


The physical nature of tokamak perturbations that would give the stochasticity 
displayed by the simple map is unclear. However, we should like the 
perturbations to be more physical in their characteristics. This can be done by 
adding terms to the map that represent low-M, low-N perturbations as well as 
terms to that would represent a diffuser coil. 


Low-M, low-N perturbations 


Such perturbations arise naturally. An area-preserving map can be defined as 


oA oA 
©,  Tg—k 2 — Yı Ye tk 
1 


with A= A(z,, yo). 
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A(x,y) is called the generating function for the map. If A,(x,y) is the 
unperturbed generating function then 


Ale, y) = Ay +48(a?—y*) cos Ge). (13) 
N, 


with J the iteration number. This perturbation adds the effect of two 
quadrupoles with opposite helicity of toroidal mode number +1. 


High-M, high-N perturbations due to localized coils 


The derivation by Martin & Taylor (1984) of the high-M and high-N 
perturbation due to localized coils gives 


B,—be “sing, (14) 
B, = b e™” cosky. (15) 


Here x = M/r, where M is the poloidal mode number of the perturbation and r 
is the minor radius of the tokamak. These perturbed fields are practical only in 
the limit of high toroidal mode number N for a finite-aspect-ratio tokamak 
(Martin & Taylor 1984). Thus, by choosing sufficiently large x, we can represent 
high-M, high-N perturbations. 
The field lines are given by 
dx B, 


dy _&, 
dz B, a 
Let h(x, y) = 2 çe sin Kx 
Y Ei KB > 
dı oh 
then a E7 = kh, (18) 
dy oh 
ae By" (19) 


Since h is independent of z, it is conserved. Now assume that the coil has a z 
length / that is very small compared with the tokamak circumference; then the 
jump in a field line is 


Ly = Uy tKlh(Xy, Yo) = Lo +Ô eosin kx, (20) 
, bi 
with öz B? h(x,, Yı) hg, Yo)- 
0 
1 sin Kx 
Theref = — 1). 21 
ou Yi = Yor ma (= =) (21) 


This map will be well behaved if «|x, —,| < 1. We apply the diffuser map, (20) 

and (21), after each of the N, iterations of the regular map required to represent 

the one full toroidal circuit of the tokamak if there is only a single diffuser coil. 
4 PLA 52 
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4.4. Deposition pattern of the field lines 


There are two ways in which we can study the deposition pattern of field lines 
on the divertor plates. In the first case we start the field lines just outside the 
last closed surface, where they undergo many circuits of the tokamak. In the 
second case we start the field lines inside the last closed surface, and radially 
diffuse them outwards after each iteration. We illustrate both methods when we 
discuss the application of the simple map. 


5. General theory of maps 


We have used maps to represent the features of magnetics of divertor 
tokamaks. Here we give a brief discussion of some relevant aspects of the 
general theory of maps. 

Magnetic field lines are the trajectories of a single-degree-of-freedom 
Hamiltonian A(x, y, £). An arbitrary trajectory x = x(2, Yo č), Y = Y(Xo: Yo, E) is 
a Hamilton trajectory, 


dx 0H 


dy OH 
de OE? (23) 
if and only if it is area-preserving, 
A(x, y) 
= 1, 24 
Öz, Yo) 
where the Jacobian 
Oxy) _ Ox Oy Ou dy (25) 


A(X, Yo) > Öz, OY URIN 
The trajectories change radically if the Hamiltonian structure (magnetic flux 


conservation) does not hold. 


5.1. Generating functions 
A map is area-preserving, 0(2,, yı)/0(24, Yo) = 1, if and only if 


aA, 
Ly = Ly İY, (26) 
oA 
and =y,+k—, (27) 
Yı Yo Öz; 


where A,(2,, Yo) is the generating function for the map. The connection between 
Hamiltonian mechanics and the topology of the map is determined by the map 
parameter k. In the limit of vanishing k, A, becomes the Hamiltonian. If A, is 
independent of iteration number, it is conserved by the map: 


Ay, Y1) = Aly, Yo) OR). (28) 
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For small k one can pick A(z, y) so that either the constant-A surfaces have 
the desired shape and/or the Hamiltonian A(x, y) has the desired shear. 


5.2. Explicit, reversed and continuous maps 


A map is explicit if the equation 
04/2, Yo) 
i 0 by, 


can be solved analytically, x, —2x,(2,,y,). Implicit maps require more 
computation. The simple map is an example of the explicit map. It is possible 
to reverse any map: 


OA(2,, 

m= yı ke, (29) 
1 
OA 


It is important to remember that for the backwards map to be explicit an 
analytical solution yọ = y,(2,, yı) must exist. In our investigation we intend to 
use only explicit maps. 

As mentioned earlier, we want to be able to determine the exact toroidal 
location of the field line or plasma particle. This requires a continuous analogue 
of the map. Any area-preserving map can be made continuous by letting 


(8) = me (31) 


WE) = yet HES, (32) 


with O<£ <1. The plane determined by the discrete map or =0 or 1 
is called the principal plane. x, and y, must lie in the principal plane. Let 
U(E) = %(%, Yo, É) and y(E) = y(%, Yo, E); then x,  x(6,) + x(x, Y1: 2—61). The 
continuous map is equivalent to a Hamiltonian that has a localized perturbation 
in each principal plane. 


6. Monte Carlo calculations of the particle motion 

Our aim is to determine the patterns of heat deposition on the divertor plates 
of a tokamak. To achieve this, we need to calculate the trajectories of the 
plasma particles in a divertor tokamak with a stochastic scrape-off layer. 
Plasma particles enter the stochastic scrape-off layer from the main plasma, 
where the magnetic surfaces are closed. Collisional effects cause the particles to 
diffuse across magnetic surfaces. The relevant kinetic equation is the Boltzmann 
equation with Lorentz collision operator: 


FF. = (34 Fa) FOUN) (33) 


4-2 


ma “re 
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Here l is the distance along the magnetic field line, C(f) is the collision operator, 
v is the velocity parallel to B, D is the spatial diffusion coefficient, 7 is 
temperature, x and y are the rectangular co-ordinates in a plane normal to 
toroidal direction, and 


T 
Cif) = blow (34) 


is the collision operator. Here u is a flow velocity. Now the collision and 
diffusion operators can be represented by their Monte Carlo equivalents (Boozer 
& Petravic 1981; Punjabi et al. 1990). After a time step of length 7, the particle 
velocity is changed by collision from v, to v,, with 


Vy = Yy— VT(v,— uu) + 6 n). (35) 
m 
The position is changed by collisions from (24, yo) to (%,,y,), with 
a, = %+(Dr), (36) 
Yı = Y+ (Dr). (37) 


The Monte Carlo procedure starts a particle with velocity v, in one of the 
principal planes of the map. After a time step 7, the particle has moved a 
distance | = vo T, so the field line map is iterated //L times, with L = 27R/N, 
being the distance associated with one iteration. N, is the number of iterations 
of the simple map that are equivalent to one circuit of the tokamak. The 
fractional part of an iteration is the £ that was introduced to make the map 
continuous. The collision and diffusion Monte Carlo operators are applied, so 
one has the new values of x, y and v and one knows the fractional distance £. 
The particle located at (x, y, £) moves a distance / = vr along the line in the next 
time step. We can use the backwards map to find the location of the field line 
in the previous Poincaré plane, and then iterate the map a distance //L +. 
Thus we can follow the particles from the plasma interior to the divertor plate. 
The electric sheath at the divertor plate implies that the electron with velocity 
v should be reversed at the plate unless mv? > ed,. The energy of the particle 
when it strikes the divertor plate can be calculated directly. To implement 
energy scattering as well as evolution of an arbitrary initial distribution 
function, the method developed by Boozer & Petravic (1981) and Punjabi and 
Boozer (see Punjabi ef al. 1990) can be used. 


7. An application of the simple map: stochastic broadening of the 
scrape-off layer 


We have used the simple map to study stochastic broadening of the scrape- 
off layer in a single-null divertor tokamak. Stochastic broadening is calculated 
from the width of the footprint of particle trajectories on the divertor plates. 
We locate divertor plates in the plane y = 1. Plasma particles following the field 
lines are led into the stochastic region by diffusion, where they strike the 
divertor plate. The footprint on the divertor plate consists of one or more 
curved stripes (Punjabi et al. 1992). The width of the stripe on the divertor 
plates is uniquely determined by the stripe formed by field lines as they cross 


Tokamak divertor maps 101 


ns. oe 
Shee eet ay'ı şa mer oT 


-0-4 0 0-4 


Fiaure 2. Application of radial expansion when k = 0-6 and 6 = 10-5. (a) represents 2048 
trajectories passing through x = 0, y = 0-1 in the principal plane. (b) and (c) represent (a) in 
the principal plane after 3-5 x 109 and 7 x 10° applications of the radial expansion operator. 
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Ficure 3. Footprint of trajetories on the divertor plate. Subscript s stands for ‘strike’. 
(p,,%,) are the g and x co-ordinates of the strike points on the divertor plate. Trajectories 
start on the surface passing through x = 0, y = 0:5, with k = 0°6 and ô = 10-11. The footprint 
is a single stripe. The divertor plate is located in the plane y = 1. 


the plane y = 1, which passes through the X-point. The toroidal location ¢ at 
which field lines cross the (x,¢) plane is determined using the continuous 
analogue of the discrete map, (6) and (7) above. The behaviour of field line 
trajectories that cross the last confining surface at a random point can be 
assessed by adding a small area-expanding step to each iteration of the discrete 
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FIGURE 4. Footprint of trajectories starting on the surface going through x = 0, y = 0°5 
for k = 06 and 6 = 10-5. Now the footprint is made of helical stripes. 
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FIGURE 5. A typical Poincaré surface-of-section in the principal plane. The trajectory starts 
at x = 0, y = 0-5, with k = 0-6 and 6 = 10. This figure represents 5 x 10° iterations. 


map, aew = (1 +Rô) yg, with R a random number between 0 and 1, and ô the 
expansion coefficient. For k = 0:6 we find that for ô < 10-3 the strike points 
form helical stripes that are independent of ô. For larger values of ô the strike 
points form a single stripe, which has a width that scales as i. From the 
physical point of view, we want this stripe to have as large an area as possible, 
so that the heat flux density is a minimum. This means that we should like the 
magnetic flux from the stochastic scrape-off layer striking the divertor plates to 
be as large as possible. 
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Ficure 6. (a) Footprint of the strike points on the divertor plate when 8 = 0. The trajectories 
start in the scrape-off layer at x, = 0 and 0:977 < y, < 1. (b) When the footprint for ô = 0 is 
unfolded in the ¢ direction, we get a sharp, narrow helical stripe. 


7.1. Stochastic broadening 


To simulate plasma-particle trajectories in a single-null divertor tokamak, we 
start with 2048 particles on the flux surface passing through (x = 0, y = 3). The 
trajectories of these particles are evolved using the Hamiltonian equations. At 
the end of each toroidal circuit, the radial expansion operator is applied. Figure 
2 shows the application of the radial expansion operator. The (x, ø) co-ordinates 
of the intersections of trajectories with the plate are recorded until half of the 
particles have struck the plate. This gives the footprint of the stochastic layer 
on the plate. The width of this footprint is defined as the standard deviation in 
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Figure 7. Scaling of average width <w) with radial expansion coefficient 6. The error in 
the mean is less than 1%. 
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FIGURE 8. Plot of ¢, as a function of Ay = |y — yol. The trajectories start just outside the last 
good surface through x = 0, on the far side of the X-point. Here k = 0-6 and ô = 0. The solid 
line represents the ¢ co-ordinate of the strike point, ¢,, as a function of Ay. The dashed line 
is the fit g, = 0-24A-928, 


the x co-ordinate of the strike points. Stochastic broadening is the width in the 
limit of vanishing ô. We find that when ô > 10-*, the footprint is a single stripe. 
For 6 < 10-5, the strike points form helical stripes. Figure 3 shows the single 
stripe, and figure 4 shows helical stripes. Figure 5 shows a typical Poincaré 
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FIGURE 9. (a) The points in the helical strip corresponding to numbers of poloidal circuits, 
N,, between 1 and 2 exhibit a concentric ring-like structure. We reverse-integrate a very 
large number of trajectories, starting on the divertor plate on the positive-x side. The 
trajectories start on uniformly distributed grid points in a rectangle exactly accommodating 
the helical stripe. This plot shows the grid points whose reverse trajectories strike the 
divertor plate on the negative-x side with 1<N, <2. (b) Grid points whose reverse 
trajectories strike the other side of the plate with 2 < N, <3. 


2; 


surface-of-section for one trajectory. To obtain the footprint in the limit as 
6-0, we put 6=0 and start trajectories in the scrape-off layer near the 
X-point. When unfolded in the ¢ direction, the footprint is a narrow and sharp 
helical stripe (figure 6). For the single-stripe case the average width of the 
footprint varies as ô: 

Kw) o ot. (38) 
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Ficure 10. Graph of n versus N,. We integrate a very large number of trajectories 
backwards. They start on the divertor plate on the positive-x side on uniformly distributed 
grid points in a rectangle that exactly accommodates the helical stripe. Solid squares 
represent n as a function of N,, where n is the number of trajectories that strike the other 
side of the plate after a fixed number N, of poloidal transits. The solid line is the fit, 
n = 37221N5'8. 


This is shown in figure 7. For ô < 107"? the average width is independent of ô. 
The & scaling is the result of the scaling Ay oc g-i (figure 8). Field lines in the 
stochastic layer spend the majority of time in the proximity of the X-point. 

We generate a rectangular mesh on one side of the divertor plate. We consider 
the nodes of this mesh as the strike points of trajectories on the divertor plate. 
Now we use the backward map, (9) and (10), and its continuous analogue to 
calculate the starting points of these trajectories on the other side of the 
divertor plate. We also record the number of poloidal circuits it takes each 
trajectory to traverse from one side of the divertor plate to the other. From this 
exercise, we observe a number of interesting features. The helical stripe is 
actually the footprint of the trajectories that circle the magnetic axis at least 
once. Within this stripe, the points corresponding to different values of poloidal 
transits show concentric rings. We denote the number of poloidal transits by 
N,- In figure 9 we show the rings corresponding to N, = 2 and 3. In figure 10 we 
plot the number n of open field lines, for a given fixed number N, of poloidal 
transits. From this plot, we see that 


no Ny. (39) 


From this, we conclude that the area of the helical stripe varies as the inverse 
square of the number of poloidal transits for the open field lines. 

In figures 11-13 we present some data on the distribution of strike points 
within the helical stripe. In the x direction the distribution is skewed towards the 
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Figure 11. (a) Distributions of strike points in the helical stripe: (å) in the x direction; 
(b) in the g direction; (c) transverse to the stripe. 
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Figure 12. For caption see facing page. 
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FIGURE 13. Dwelling time in the poloidal direction 6, as a function of Ay. 


X-point (figure 11a). In the ¢ direction the distribution is skewed towards 
increasing ¢ (figure 11b). Transverse to the stripe, the distribution is normal 
(figure 11c). We also observe that the dwelling time in the toroidal direction for 
the open field lines before they strike the divertor plate has a self-similar 
pattern as shown in figure 12, while the dwelling time in the poloidal direction 
has a discrete structure as shown in figure 13. 


7.2. Summary of results 


Particles that diffuse across the last confining surface of a tokamak with a 
divertor can follow open magnetic field lines to the divertor plates. The strike 
points of these field lines on the divertor plates lie in discrete stripes. The region 
between the last confining magnetic surface and the separatrix of the ideal 
divertor is called the stochastic layer. The open magnetic field lines of the 
stochastic layer enter the region from one divertor plate and leave, striking 
another divertor plate, after an integer number N, of complete poloidal circuits. 
The discrete stripes on the divertor plates correspond to the number of poloidal 
circuits made by open field lines. The area of the stripes scales as N;”. The open 
field lines that approach the last confining surface make many poloidal circuits 
before striking a divertor plate in a stripe that is much narrower than the 
stochastic layer ( 15%) but none the less of non-zero width. An open field line 
that makes many poloidal circuits in going from one divertor plate to another 
in essence densely fills the stochastic layer except for the regions in the latter 
that are occupied by magnetic islands. 


Frcure 12. Dwelling time for trajectories in the toroidal direction ¢, as a function of distance 
Ay = |y—y,| from the last good surface. The trajectories start on the line x = O in the scrape- 
off layer. We determine the number of toroidal circuits. N,, that it takes the trajectories to 
reach the divertor plate as a function of distance Ay. (a) Ay = 1074; (b) 1077; (e) 10-9, 
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7.3. Discussion of results 


In an axisymmetric divertor the spreading of the heat is due to plasma 
transport, which is often modelled by a diffusion coefficient D of order 1 m? s7?. 
During the time 7 that the plasma takes to flow to the divertor chamber, the 
plasma diffuses across open magnetic field lines that contain a flux proportional 
to (Dr). The open field lines that contain flowing plasma have the same 
magnetic flux near the main body of the plasma and near the divertor plates. 

By adding asymmetric magnetic perturbations, one can widen the region of 
heat deposition on the divertor plates. However, we find that the magnetic flux 
of the open field lines in the stochastic region is a factor of roughly six larger 
than the width in flux of these field lines when they strike the divertor plate. In 
effect, the open field lines of the stochastic region encircle the plasma roughly 
six times before striking the divertor plates. The relatively thick region of 
plasma flowing along open field lines in the stochastic region may provide an 
important shielding of the main plasma from impurities. 

Magnetic field line trajectories are the trajectories of a 1}-degree-of-freedom 
Hamiltonian (Boozer & Rechester 1978; Cary & Littlejohn 1983). Consequently, 
the tokamak separatrix and divertor are archetypal for the behaviour of a 
region of bounded Hamiltonian trajectories surrounded by a region of open 
trajectories. Plasmas flow rapidly along magnetic field lines and diffuse slowly 
across them. The study of magnetic field line behaviour and the closely related 
plasma behaviour are important not only for their application to tokamaks but 
also for their application to other Hamiltonian, or near-Hamiltonian, systems. 

The behaviour of field lines near a tokamak separatrix has been studied 
previously. Early work (Boozer & Rechester 1978) used the Chirikov island 
overlap criterion, which is a reasonable approximation for calculating the 
stochastic region but overestimates the spreading of heat on the divertor plates. 
Recently, Pomphrey & Reiman (1992) integrated magnetic trajectories in a 
field produced by wires. They noted that the extremely thin spiral regions of 
impact of magnetic field lines started near the last closed surface on the divertor 
plates. 

The major advantage of our method is that field line trajectories can be 
followed many orders of magnitude faster with a map than with an explicit 
integration. The speed of map iterations allows one to study the structures 
formed by the perturbations in more detail and to investigate the effects of 
qualitatively different types of perturbations. 

The complicated filling of the stochastic region by the open field lines gives 
an associated complexity to the ambipolar electric potential that exists along 
any open magnetic field line that is embedded in a plasma. Although we have 
not calculated the particle transport from the E x B drifts associated with the 
ambipolar potential, we expect a diffusive effect that would smear the stripe 
structures. 


This work is supported by United States Department of Energy Grants DE- 
FG05-88ER53265 and DE-FG05-90ER54106. 
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